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Abstract. We study the emergence and dynamics of pointer states in the 
motion of a quantum test particle affected by collisional decoherence. These 
environmentally distinguished states are shown to be exponentially localized 
solitonic wave functions which evolve according to the classical equations of 
motion. We explain their formation using the orthogonal unraveling of the master 
equation, and we demonstrate that the statistical weights of the arising mixture 
are given by projections of the initial state onto the pointer basis. 



PACS numbers: 03.65.Yz, 47.45.Ab, 05.40.Jc 
1. Introduction 

The influence of environmental degrees of freedom has been identified as the key 
concept in explaining the classical behavior of macroscopic systems in a quantum 
framework [1-3]. According to this notion a preferred set of localized system states - 
called the pointer basis [4] - is induced in the course of the interaction of the system 
with its surrounding. Most characteristically, any initial superposition of these pointer 
states gets rapidly mixed, while the only states who retain their purity for a long 
time are the pointer states themselves. While the basic ideas behind the decoherence 
process seems to be settled, it still remains an open problem to understand the 
emergence, the dynamics, and the main properties of the pointer states for microscopic 
realistic environments. 

Several strategies have been proposed so far for determining the pointer basis 
given the environmental coupling. In [5] the suggestion was made to sort all pure states 
in the Hilbert space according to their linear entropy production rate. The pointer 
states are then identified with the states having minimal loss of purity. Similar results 
are obtained by the approach of [6-8] which is based on a time evolution equation 
whose solitonic solutions are identified with the pointer states. So far, this concept 
has been applied to the damped harmonic oscillator [7, 8] and to a free quantum 
particle coupled linearly to a bath of harmonic oscillators [6,8]. There, the solitonic 
solutions of the corresponding nonlinear equation are coherent states and Gaussian 
wave packets, respectively. Moreover, the decoherence to Gaussian pointer states was 
proved to be generic for linear coupling models [9] . 

In this paper, we extend the analysis from linear models to a non-perturbative 
treatment of the interaction. We focus on the model of collisional decoherence, which 
provides a realistic description of the decoherence process generated by an ideal gas 
environment. Notably, experiments with interfering fullcrene molecules displayed a 
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reduction of interference visibility in agreement with this model [10-14]. We derive the 
corresponding pointer states which are shown to form an ovcrcomplete, exponentially 
localized set of basis states, and we prove decoherence to these states using the 
orthogonal unraveling of the master equation [15,16]. This stochastic process on the 
one hand provides the statistical weights of the pointer basis, and on the other hand 
represents an efficient way of solving master equations which exhibit pointer states. 
Moreover, it explains the emergence of classical Hamiltonian dynamics. 

The main result was already announced in [17]. Here, we provide a more detailed 
explanation of the proofs and derivations. As an extension of [17], we prove the 
decoherence dynamics for a more general situation where the localization rate of 
collisional decoherence is in a non-saturated regime, and we utilize the relative entropy 
in order to illustrate the emergence of the statistical weights. 

While the above results are derived within the framework of decoherence theory, 
they can also be applied to dynamic reduction models which propose a modification 
of the Schrodingcr equation by means of nonlinear and stochastic terms. In fact, 
the observational consequences of the Ghirardi-Rimini- Weber (GRW) spontaneous 
localization model [18, 19] are equivalent to the ones of collisional decoherence, since 
they are described by the same master equation [20]. The present work therefore 
applies to the GRW model. In particular, it provides the corresponding pointer basis. 

The structure of the article is as follows. In Sect. 2, we briefly review the notion 
of pointer states, and we summarize the method for determining the pointer states 
discussed in [6-8]. In order to motivate the approach, we consider a two level system 
and a dephasing process. This method is then applied to collisional decoherence 
in Sect. 3, which provides a set of solitonic states to be regarded as 'candidate' 
pointer states. We show in Sect. 4 that these solitons form an ovcrcomplete basis 
of exponentially localized states, and we give an expression for their spatial extension. 
Moreover, we demonstrate that the 'candidate' pointer states move on classical phase 
space trajectories if they are sufficiently localized. In Sect. 5, we briefly review the 
method of quantum trajectories, focusing in particular on the orthogonal unraveling of 
the master equation. This stochastic process is then applied to collisional decoherence 
which allows us to show that the 'candidate' states are indeed pointer states in the 
sense of the definition given in Sect. 2. Furthermore, we use the orthogonal unraveling 
to show that the statistical weights of the pointer states are given by the overlap with 
the initial state. We present our conclusion in Section 6. 

2. The pointer basis 

2.1. Definition of pointer states 

To motivate the definition of pointer states, let us consider the quantum dynamics of 
the damped harmonic oscillator. Its evolution can be described by a master equation in 
Lindblad form defined by the standard Hamiltonian H = fiwa'a, and a single Lindblad 
operator L — a, with associated rate 7 [21]. We take as initial state a superposition of 
two quasi-orthogonal coherent states 

[V'o) = ci|a ) + c 2 \f3 ) , with |a - f3 \ 2 > 1 , (1) 

which satisfy a|a) = a\a), with aeC. It is then easy to show that for times larger 
than the decoherence time tdoc = 2 \a — /3q| -2 /"f the solution of the master equation 
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is well approximated by [21,22] 



pta | Cl | 2 |a t )(a t | + |c 2 | 2 |/3 t )(A| 



if t > t d , 



(2) 



with a t — ao exp [—iwt — r yt/2). Thus, any coherent state remains pure during the 
damped time evolution, while any superposition of distinct coherent states decays 
into a mixture (with a decay rate 7d oc = lAdec ^ 7) whose statistical weights are 
determined by the initial overlaps KaolV'o)! and | {/?o I V^o) I ■ Due to this property, the 
coherent states are to be identified with the pointer states of the damped harmonic 
oscillator. 

The above observation serves as the starting point for the following definition 
of the pointer states for an open quantum system evolving according to a Lindblad 
master equation dtp = Cp. One says that the system exhibits a pointer basis if its 
dynamics involves a separation of time scales, characterized by a fast decoherence 
time tdec, such that for any time much larger than t& cc , the evolved state is well 
approximated by a mixture of uniquely defined pure states P Q = |7r Q )(7r Q | which are 
independent of the initial state poj 



with Prob(a|po) 0, J da Prob (ot\po) = 1. Following the above example, we further 
demand that for initial states po, which are superpositions of mutually orthogonal 
pointer states P^, the probability distribution Prob (ot\po) = Y]g wpS (a — (3) is given 
by the initial projections 



The pointer states P Q initially form an (overcomplete) basis, and they may evolve in 
time, though slowly compared to t& cc . 

The name pointer state was coined in [4] due to its relevance for the physical 
description of a measurement apparatus. A measurement device which probes an 
observable A is constructed such that macroscopically distinct positions of the pointer 
or indicator are obtained for the different eigenstates of A. For a quantum system 
initially prepared in an eigenstate of A the read-out will display the corresponding 
eigenvalue with certainty provided these pointer states remain pure during the time 
evolution. On the other hand, if the quantum system is prepared in a superposition 
of eigenstates of A, we expect the pointer not to end up in a superposition of different 
read-out states, but rather to be at a definite position, though probabilistically, with 
probabilities given by the overlap ((4|). 

We emphasize that the importance of pointer states goes beyond the physics of 
measurement devices and the quantum-to-classical transition since they are also a 
practical tool for the solution of master equations. Knowing the pointer states P a , 
their time evolution P a (t), and their probabilities Prob(a|/9o), one can immediately 
specify the solution of the master equation for any initial state and times greater than 
the decoherence time. Since the decoherence time is generically much shorter than 
the system and dissipation time scales of the pointer state motion, this allows one to 
capture a large part of the system evolution without solving the master equation. 

2.2. Pointer states of pure dephasing 

A practical way to obtain the pointer states P Q , for a given environmental coupling, 
was discussed in [6-8] . We will illustrate this method by means of a two level system, 





(3) 



w /3 = Tr(p P/3(0)) . 



(4) 
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suspect to a dephasing environment. The corresponding master equation in interaction 
picture, 

dtP=l{o- z pa z - p) , (5) 

is characterized by the Lindblad jump operator L =<t z ^/j (7 > 0). In Bloch 
representation, pt = 1/2 (I + a (t) ■ er) with Bloch vector a (t) — Tr (crp (t)) and Pauli 
matrices er = (a x , o~ yi o~ z ), the solution reads 

a (t) = (e- 2 ^a x (0) , e" 2 ^ (0) , a z (0)) . (6) 

Thus, the Bloch sphere {a : \a\ — 1}, which represents the set of all pure states, is 
projected onto the z-axis in the course of the dephasing process. This implies that the 
decohered state a (00) = (0, 0, a z (0)) is a mixture of the eigenstates of a z (denoted 
by Pi = I 4X1 I and P T = I TXT I respectively), 

Poo = Tr[P TA) ] P T +Tr[P iPo ] P x . (7) 

The comparison with ([3]) shows that the north and the south pole of the Bloch sphere 
(corresponding to Pj_ and Pf ) form the pointer basis of the pure dephasing process. 

Since the solution of the master equation will not be at hand in general, one 
requires a method which yields the pointer states without the knowledge of p t ■ To 
motivate this, we notice that the north pole of the Bloch sphere is the asymptotic 
end point of the trajectory illustrated by the thick line in Fig. 1. This trajectory is 
generated by an equation of motion with the following properties, (i) It is nonlinear 
because it must distinguish pointer states from their superpositions, (ii) It preserves 
the purity of pure initial states, i.e. an initial state which lies on the Bloch sphere 
remains on the surface. (Hi) The generated trajectory follows the exact solution of 
the master equation as close as possible. In order to find such an equation of motion, 
it is suggestive to minimize the distance of the initial increments, i.e. 

min \C (a) — dta\ , with a ■ dta = . (8) 

dta 

Here, C denotes the generator of the master equation ([5]) in Bloch representation and 
d t a is subject to the condition that the generated trajectory remains on the surface of 
the Bloch sphere. The solution of the above optimization problem reads, in spherical 
coordinates, 

(f,<p,0) = (O,O,- 7 sin(20)) . (9) 

Since the sine is positive for 9 £ (0,7r/2), the solutions of these equations tend 
asymptotically towards a pointer state of the system, see Fig. 1. The equator of 
the Bloch sphere forms a set of unstable fixed points of ^ . 



2. 3. Nonlinear equation for pointer states 

Let us now generalize the above argument to general Markovian master equations. 
Replacing the Euclidean norm by the Hilbert-Schmidt norm || • ||hs m the space of 
operators, the generalization of ([5]) to higher dimensional systems reads 

min\\C(P)-d t P\\ 2 m , (10) 

dtP 

where the minimization is with respect to all evolution equations dtP = f (P) which 
propagate P within the set of pure states, such that P 2 = P* . It can be shown 




Figure 1. Bloch representation of a two level system subject to pure dephasing, 
as described by J5}. As t —* oo, the initial state po = P is projected onto the 
z-axis, implying that the poles (Pj and Pj ) are the pointer states. The thick 
line indicates a trajectory within the set of pure states which connects the initial 
state to a nearby pointer state (the north pole). The equation of motion for this 
trajectory has minimal distance from £(P) among all equations which generate 
pure state trajectories. 



that the general structure of an equation that evolves state vectors preserving 
their normalization has the structure d t \ip) = (A,/, — {ip\Aj,\ip) + B,/,) \tj)}, with in- 
dependent, hermitian and anti-hermitian mappings = A^ and B^, = — B^ . This 
implies that the equation of motion for the projector P = |f/')(V ; l must be of the form 
dtP = [P, [P, Xp]] , where Xp := Ap + [Bp, P]. Using this form, the optimization problem 
(|10P reduces to minx \\£ (P) — [P, [Pj X]] || HS . As shown in [7] the solution is determined 
by the generator of the master equation, X nun = C(P). Hence, the generalization of 
© reads as [6-8] 

dtP= [P,[P,£(P)]]. (11) 

Motivated by the example in Sect. 2.2, one conjectures that the asymptotic solutions 
of pip provide the pointer states in more complex systems as well. 

It will be important below that Eq. is known also in another context: it 
corresponds to the deterministic part of the orthogonal unraveling of the master 
equation. As we will demonstrate in Section 5, one can use this specific unraveling 
to prove for a specific model that the asymptotic solutions of (fTTj) indeed provide the 
pointer states. 

3. Pointer states of collisional decoherence 

3.1. Collisional decoherence 

In order to assess the nonlinear equation in the context of a nontrivial 

environmental coupling we now apply it to the model of collisional decoherence 
[23,24]. The latter describes the motion of a quantum test particle in an ideal gas 
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environment and it accounts for the quantum effects of the scattering dynamics in a 
non-perturbative fashion. The corresponding master equation has Lindblad form, 

d t p = - ~ [H,p] + Jdq (l q pC q - ±C q L qP - \ P ^)j , (12) 

where the jump operators are proportional to momentum kick operators, L q = 
\AyG (q)e lqx (with position operator x). The continuous label q has the meaning of a 
momentum transfer experienced by the test particle with G (q) > the corresponding 
distribution, J dqG(q) = 1; 7 is the collision rate of the gas environment. The Id 
equation of motion thus reads 

d t p = - [H, p] + 7 / dqG (q) e^/ h pe-^' h - 1P . (13) 
in> J-00 

It leads to a localization in position space, i.e. to a loss of spatial coherence, as can 
be seen by switching to the interaction picture, p — e Mt / h p e ~ Mt / h ^ anc j the position 
representation, 

d t (x\p\x') = — F (x — x') (x\p\x'} . (14) 

The decay rate of the spatial coherences is thus characterized by a localization rate 
F (s) > which is related to the momentum transfer distribution G (q) by 

F(s) =7(1- r dqG(q)e^ h ) . (15) 



Since the Fourier transform of the distribution G (q) tends to zero for large distances s, 
the localization rate saturates for large s at the maximum value given by the collision 
rate, F (s — > 00) = 7 , which can be interpreted as the limit where one collision is 
sufficient to reveal the particles 'which path' information. This behavior is in sharp 
contrast to linear models, where the localization rate grows quadratically, and thus 
approaches infinity in the limit of a large separation s. 

3.2. Determining the pointer states of collisional decoherence 

In order to apply the nonlinear equation (fTTjl to collisional decoherence fp~3|) of a 

free particle, that is H = p 2 /2m, we rewrite the projector equation (llip in vector 
representation and choose for C the Lindblad form (fT2")l . which gives [7,8] 

W) = ^(H-(H)^)|V> 

+ J°° dq (Vj).0 (L q - (L,)*) - \ (LjL, - (L+L fc >*)) |V) • (16) 

The expectation value (H)^ is disregarded in the following, since it contributes only 
an additional phase. Now, we choose the jump operator of collisional decoherence, 
Lq = \J 7G (q)e lqx , and switch to position representation, which yields 

dt^t (x) = - ^-d 2 x A (x) + in (?) A [iV-tll (x) , (17) 
Ami L J 



A 



\M 2 ] {*) = 7 (\M 2 *G(x) — j°° dy\M 2 (V) (iVtl 2 * G) {y)j . (18) 



Here, g * h(x) = J_ dyg (y) h (x — y) denotes the convolution and G (x) is the Fourier 
transform of G (q) , i.e. G (x) = dq G (q) exp (iqx/H) . 
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Figure 2. Formation of pointer states: an initial superposition of counter- 
propagating, localized states is evolved numerically according to the nonlinear 
equation 119> . It forms into a solitonic solution which moves with a fixed shape 
and constant velocity. These solitons are interpreted as the pointer states of 
collisional decoherence. 



The two summands in (jTTJ) have counteractive effects on the temporal evolution 
of the wave function: the coherent term leads to its dispersion, whereas the 
second, incoherent summand tends to localize the solution. In order to explain this 
localization, we note that the second summand in (|18p is independent of x. This 
implies that the centered parts of the wave function, where the convolution |V't| 2 *G' (x) 
exceeds the constant term in (TT51) . get amplified, i.e. dfipt > Oj whereas the tails of the 
wave function get damped, i.e. dttyt < 0. As a consequence of these competing effects, 
solutions of (fT7j) evolve towards solitonic states where both effects are in equilibrium, 
such that the state moves with fixed shape and constant velocity. As discussed above, 
these solitons are candidates for the pointer states of collisional decoherence. 

Assuming the momentum transfer distribution G (q) to be a centered Gaussian 
with variance o~q , we can rewrite (|17|l in dimcnsionless form. 

dr^r (y) = - T d 2 y tp T (j/) + if T (y) / dy'\ip T (y r ) | 2 

x (V(^o 2 /2 _ y°° dy "^ T |2 e -o/-</') 2 /2^ (19) 

Here we use the dimensionless variables y = aox/h and r = yt to define the 
dimensionless wave function ip T (y) = y/ h/aaipr/i (%/cg) ■ Notably, Eq. (fT9|) 
depends only on the single dimensionless parameter k = Oq / (mtvy) . 

Figure 2 shows a numerical solution of (fl9|) , where we choose a superposition of 
two counter propagating localized states (f>i^ as the initial state, ipo (x) = c\<j>\ (x) + 
C2<f>2 (x). As expected from the above discussion, the (modulus of the) solution 
converges to a soliton. Moreover, we find that the soliton inherits its initial position 
and momentum expectation value from that localized component of the initial 
state which has the greatest weight Cj, |cj| > |cj - ^<|. Similar observations are found for 
various other initial states. 
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4. Properties of the soliton basis 

We proceed to characterize the solitonic solutions of (fl~7|) . In Sect. 4.1, the 
consequences of the conservation of probability on the phase of the solitons are 
analyzed, allowing us to predict the asymptotic shape of the solitons in Sect. 4.2. 
In Sect. 4.3 we estimate the spatial extension of the solitons, followed by the proof 
that they form a basis of the Hilbert space in Sect. 4.4. Finally, in Sect. 4.5, we discuss 
the dynamics of the solitonic solutions in the presence of an external potential. 



4-1. Consequences of the continuity equation 

As demonstrated in the previous section, the nonlinear equation (|17[) exhibits solitonic 
solutions 7r t (x) in the sense that the modulus of n t (x) moves with constants shape 
and velocity, i.e. 

7r t (x) = f(x-vt)e^ t \ (20) 

with / and g real. In this section, we analyze the general structure of the phase 
g (x, t), which will be relevant subsequently. The time derivative of a solution \tp t (x)\ 2 
of ([T7|) . yields the continuity equation for ip t (x), 

d t \^ t (x)\ 2 = --d x lm(r t d x ^t) + 2\^ t (x)\ 2 A[\^ t \ 2 ] (x) . (21) 
m 

Plugging the solitonic solution (|20p into (|2"Tj) , gives 

- 2A [f 2 ] (x - vt) - vd x log f 2 (x-vt) = - - {dig (x, t) 

m 

+d x g(x,t)d x \ogf 2 (x-vt)) . (22) 

Here we have used that A [/ t 2 ] (x) — A [f 2 ~\ (x — vt), which follows from ft (x) = 
f (x — vt) . The time dependence of the left hand side of (|22|) corresponds to a 
spatial shift. Thus, also the right hand side of (|22|) must exhibit such a simple time 
dependence, which implies that 

-vd x r(x,t) =d t r(x,t) , (23) 

where r (x, t) denotes the right hand side of (f2"2"l) . It follows that 

- vdlg (x, t) - vd 2 x g (x, t) d x log f 2 (x - vt) = d t d 2 x g (x, t) 

+ d t d x g (x,t) d x \og f 2 (x - vt) . (24) 

Since this equation must hold for all x, v and t, we may assume that the equality holds 
already for the summands, such that 

-vd 2 x g(x,t) =d t [d x9 (x,t)} . (25) 

Therefore, the temporal and spatial dependence of the phase has the general structure 

g(x,t) =<t>(x-vt)+ X (t) , (26) 

with unknown functions (f> and x- 
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4-2. Asymptotic form of the solitons 

To explore the tails of the solitonic states ir t (x) let us consider the form of (fT7|) for 
asymptotically large positions. It reads 

dtipt (x) ~ - Tr—dltfjt (x) ~ jcu/tpt (x) , for |x| — > oo , (27) 

with = dy\Tp t \ 2 (y) (lipt] 2 * Gj (y) a -0-dependent, positive constant. Inserting 
the solitonic form (|2"0j) into (|2"7|) yields 

i<9 t g (i, t)f(x- vt) = i^— d 2 x f (x - vt) - f (x - vt) (d x g (x, t)f + vf (x - vt) 

- — [dxf (x - vt) d x g (x, i)+f(x- vt) dig (x, t)] 
m 

- ia,pf (x - vt) . (28) 

Using (|26p. we find that both d x g (x, t) and d^g (x, t) are only a function of x t — x — vt, 
and accordingly, that also the left hand side of must be a function of x t . If follows 
that x (*) is & t most linear in t (that is x (*) = Xit + Xo , with unknown constants xo 
and xi). Considering the real and imaginary part of (|28[) separately, one obtains two 
coupled (second order) differential equations 

vd x f - ja^f = - (d x fd x <t> + bd 2 x A , (29) 



m \ 2 

(xi - vd x <j>) f=JL(%f-f {d x( f>} 2 ) , (30) 

where / = / (x — vt) and <f> = <f> [x — vt) . This set of equations has two unique 
solutions 

/(aO=e ±fe N, (31) 
<j>(x) = Tsg n(x)^{v+^t)x ) (32) 

where the constant k > depends on the boundary condition for (|29[) (which can be 
determined only by solving the full nonlinear equation (|17p ). The solution with the 
positive exponent in (|31|) is irrelevant, since it is not normalizable. Figure 3 confirms 
that the tails of the numerically obtained solitonic solutions of (fT9")) are in agreement 
with the functional form (|31[) ; they are straight lines in the semi- logarithmic plot. 
This shows that, unlike in linear models [25] where the pointer states are Gaussian, 
the pointer states of collisional decoherence are exponentially localized. 



4-3. Size of the solitons 

An important characteristic of the pointer states is their spatial extension. As 
explained in [17], the latter can be related to the experimentally accessible one-particle 
coherence length of a thermal gas. We will determine the pointer width in this section, 
and apply the result later, when studying the dynamics of pointer states in an external 
potential. 

As a first step, consider the standard deviation av of the numerically obtained 
dimensionless solitonic solution |7r(y)| 2 of (fT9"]l as a function the dimensionless 
parameter k = a^/ (jmK) . As shown by the solid line in Fig. 4, the size ov increases 
linearly with re . 
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Figure 3. Semi-logarithmic plot of the numerical solitonic solution of l|19|l . The 
graph clearly demonstrates that the pointer states have exponential tails. 



This observation can be reproduced by a simplified model which has the practical 
advantage that it can be applied to more complex systems, such as 3D gases with 
a microscopically realistic localization rate F [17]. The idea of the model goes as 
follows: the ideal gas environment consists of particles which collide with the system 
at a rate 7. At each collision, the ambient particles gain position information, such 
that the wave function gets spatially localized to a length scale £\ oc determined by the 
localization rate F, see ([T4|) . After the scattering event, the particle disperses freely, 
until it gets localized again by a subsequent collision. The pointer width a n is then 
obtained by averaging the time-dependent width of the wave function over the waiting 
time distribution of a Poisson process. 

More specifically, we assume that the length scale £\ oc is characterized by the free 
parameter 

F(e loc ) 



F(oo) 



(33) 



Using (|15p and taking the momentum transfer distribution G (q) to be a centered 



Gaussian with variance Oq , we obtain 



with 



dependent size 



loc 



= -21og(l- 



«lo 



«!oc) 



(34) 



The free dispersion after the collision yields the time 



0-1 (t) = 4 C + 



fit 



2m£\ oc 

Finally, the average over the waiting time distribution Prob (t) = ^e" 

(7^=1 drProb (r) - f Ata^ (t) 
Jo t Jo 



(35) 



7* 



gives 
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Figure 4. Spatial extension of the solitonic solution of l|19j l as a function of the 
dimensionless parameter k = a^/ (■ymh). The solid line represents the numerical 
solution of JT9J1. The result of the localization model (|37fl with parameter 
a loc = 0-4 is given by the dashed line. 



^ aioc V- , (36) 

where we use a linearization of o~ v (t) in the second line. The dimensionless version of 
(j36j) reads 

_ _ ovog 

0~TT — Z 

n 

= a loc + —^—K. (37) 
4a loc 

The dashed line in Figure 4 shows that the form of this equation agrees with the 
numerical solution of ([T9"l) ; the fit yields a value of ai oc ~ 0.4 for the parameter 
characterizing the localization length scale. 



4- 4- Completeness of the soliton basis 

Our next aim is to show that the solitonic solutions of (|17|) . which are interpreted 
below as the pointer states of collisional decoherence, form an overcomplete basis. 
For this purpose, we first present a general method to construct a whole manifold 
of solutions of (|11[) given a specific one. It relies on the symmetry properties of the 
corresponding master equation. Since collisional decoherence exhibits Galilean (i.e. 
translation and boost) invariance, it is then easy to show that the pointer states of 
this model form an overcomplete basis. 

Suppose there is a family of unitary operators U t , satisfying 

U t V(p) Ul =X>(UtpU;), (38) 

9 t U t = 4[H,U t ], (39) 
in 

where we denote by V the incoherent part of the master equation, C (p) = [H, p] / (ih) + 
V (p). Then, given a solution P t of the nonlinear equation (JTTJ) , also U t Pt constitutes 
a solution of (fTTj) . 
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This can be verified easily: 
[UPUt, [UPUt, C (UPUt)]] = _L (HUPU f - UPUtH) + U [P, [P,X> (P)]] U f 

= 4(HUPU f - UPU f H + UHPU f - UHPU 1 + UPHUt 

in 

- UPHUt) + U[P,[P,X>(P)]]U t 



12 



P, 



P,-[H,P]+Z>(P) 

in 



= — [H,U]PU f +U 

in 

-iuP[Ut,H] 

= d t (UPUt) t (40) 

where we dropped the time argument for brevity. Here, the first equality makes 
use of ([35]) and the unitarity of U, and the third equality is due to the relation 
[H, P] = [P, [P, [H, P]]]. In the last line, we use (JUJ and ([39]) . 

Let us now apply this to the Galilean invariance of master equation (fTB"]) . Wc 
will see that the phase space translations \J t = T s ,u = exp (i (u t x — s t p)) satisfy the 
symmetry conditions (|38[) and (|39[) provided the time dependence of St and u t has the 
particular form 



s t = uq tjm + so . 



(41) 
(42) 



The latter enact a phase space translation in accordance with the free shearing motion. 
Let us first verify condition ([38]) using T s ,uf (x) T| u = / (x — s). 



T s , u V(p)Tl 



7 



7 



/OO 
dq G (q) T s , u e^' h p e^^T^ - T^pT^ 
-OO 

/OO 
dg G (q) e^-^T^pT^ e -*«(*-)/» _ 7 T Sitt pTt >u 
-OO 



(43) 



= X>(T s , uP Tt iU ) 

In order to verify (|39[) , use the Campbell-Hausdorff formula to rewrite the translation 
operator as 



-,- f.iHX\ ( .s t p\ ( 

J s ,u = exp ^— J exp J exp ^ 

The time derivative thus yields 



. s t u t \ 
' 2h J 



(44) 



-r [ u t T s ,u 



1 



s*T Sj „p - - (u t s t + s t u t ) T SiU 



n \ m Zm 
1 

ih 



rn 

,2 



i_ T 

2m' s 



(45) 



where the shearing transformation ([4"T]) and (f4"2"]) is required in the second line. This 
confirms ([39]) for H = p 2 /2m . 

We conclude that the nonlinear equation ([T7]) exhibits a family of solitonic 
solutions Pr = T S!tl PT s u , parameterized by the phase space coordinate T = (sq,uq). 
In order to verify that this family forms an ovcrcompletc basis, let us consider a specific 
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class of phase space representations. According to [26], any Hilbert-Schmidt operator 
A can be represented as 

A= J dIM(r)T s , u QTl itl , (46) 

provided Q is a trace-class operator, i.e. < Tr (VQ^Q) < oo. Here, J dT- denotes a 
phase space integral and A (T) is a function of the phase space coordinate T. Choosing 
for A the identity I, and for Q the solitonic solution Poo of (JTTJ) with vanishing position 
and momentum expectations, we obtain a resolution of the identity in terms of the 
solitons P r = T s ,„P ,oT^ u , 

y"dIY(r)P r . (47) 

This demonstrates that the pointer states of collisional decoherence form an 
overcomplete basis. 

4-5. Dynamics in an external potential 

So far, we have characterized the solitonic solutions of the nonlinear equation (|17p 
which applies in the absence of an external force. If an additional potential is present 
the corresponding nonlinear equation contains an additional term V (x) / (ih) on the 
right hand side of (fl~7|) . The numerical treatment shows that the solutions still converge 
to localized wave packets, which, however, change their shape and velocity in the 
course of the evolution. We find that the center of these wave packets moves on the 
corresponding classical phase space trajectory provided the collision rate is sufficiently 
large. We first summarize our numerical findings and then proceed with an analytic 
explanation. 

Figure 5 shows the position and momentum expectation values of the numerical 
solution of the nonlinear equation, in case of an anharmonic external potential of the 
form V (x) — ax 4 — bx 2 , a, b > (starting from an Gaussian initial state). The panel 
on the left hand side of Fig. 5 was obtained in the limit of a vanishing collision rate 7 
(i.e. k — > 00), which turns 1|19|) into the Schrodinger equation. The solution therefore 
disperses, and the solid line shows a typical evolution of the phase space expectation 
values. The dashed line, on the other hand, gives the classical trajectory of the phase 
space point where the initial state is localized. The result for a large collision rate 7 (or 
small k) is shown on the right hand side of Fig. 5. Here, the initial state turns rapidly 
into a soliton whose expectation values move on the corresponding classical trajectory. 
This illustrates that the temporal evolution turns from quantum to classical dynamics 
with increasing collision rate 7 (i.e. decreasing k). We made similar observations with 
various other potentials. 

In order to explain the numerical observation, first consider a particle in a linear 
potential V (x) = ax. The corresponding nonlinear equation reads as 

9t*lH (») + -^—pl^ ( x ) = Jh ax ^ + ( x ) l^*! 2 * ^ ( x ) 

dy\A\ 2 (y)(\A\ 2 *G)(y)]. (48) 

As discussed in Sect. 4.1, the field- free version of this equation (a = 0) exhibits 
uniformly moving solitonic solutions of the form 

i/H (x) =f(x- vt) exp (* [<f> (x - vt) + X (*)]) ■ (49) 




Figure 5. Time evolution of pointer states in an anharmonic fourth order 
potential (solid line). The dashed line shows the corresponding classical phase 
space trajectory, (a) The collision rate 7 vanishes leading to dispersive quantum 
dynamics, (b) The collision rate 7 is large, such that the dynamics of the pointer 
state is indistinguishable from the classical trajectory. 



This implies that (|4"5|) has solitonic solutions of the form 

tpt(x) = f(x-x t )exp(ig{x-x t ,t)), 



a 



i(x,t) = 4>{x) +x' (t) - -rtx 



(50) 
(51) 



with x t = vt — at 2 /2m and x' (*) = X (*) — 2a J q drx T /h. In order to verify this 
statement, we evaluate the left hand side of (|48|) with the ansatz given by ([50]) . This 
yields 



dt^t (x) + — :d 2 ipt (x) = < 
zmi 



Gt 2d 

"' 1 jjff + ifd t x' (t) + —ix t f - v (d x f + ifd x 



h 
m 



d x fd x 4> - l -d 2 j + y (d x 



1 



(52) 



with g = g (x — Xt,t), f = f (x — xt) and 4> = (f>(x — Xt). The expression can be 
further simplified by noting that the free soliton (|49|) is a solution of the field-free 
version of implying that 



7/ A [f 2 ] (x) = ifdtx (*) -v(d x f + ifd x 



-[d x fd x <f>--d x f + -f(d x 



+ \fdU) .(•»>!) 



2 XJ T 

([5^)1 and the above definition of x' (0; one 



with A [f 2 ] (x) defined in (T5J|. Using 
finds that 

bit (x) + (x) 

which confirms that ipt (x) evolves according to (I48| . 

We conclude that in a linear potential the pointer states have the same shape 
as in the field-free case and they are uniformly accelerated like a classical particle. 



1 r 91 

— axi; t (x)+iipt(x)A ipt\ (x) , (54) 
in L J 
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For general potentials, this implies that the pointer states follow the corresponding 
classical motion, provided the spatial width of the solitons is sufficiently small, such 
that the linearization of the potential is justified over their spatial extension. Since 
the size of the pointer states decreases with the collision rate (see Section 4.3), the 
pointer states must exhibit classical dynamics in the limit of large collision rates. 



5. Orthogonal Unraveling 

Thus far, we have calculated 'candidate' pointer states as the solitonic solutions of 
(|17p . and we have studied their properties and dynamics. Next, it will be shown that 
these 'candidates' are genuine pointer states P a in the sense of (1). Moreover, we find 
that the statistical weights Prob (cx\po) of the pointer states are given by the overlap 
of the initial state po with the initial pointer state P a (0), i.e. 

Prob(aM = Tr [p P a (0)] . (55) 

In order to verify the above conjectures, let us now make use of the formalism of 
quantum trajectories [15,21,27] to solve the master equation (fl3|) . More precisely, a 
specific quantum trajectory method, the orthogonal unraveling , is distinguished by the 
physics of pointer states because the deterministic part of the associated stochastic 
differential equation coincides with the nonlinear equation (|11[) . 

We start with a general description of the quantum trajectory approach and 
the orthogonal unraveling in Section 5.1. The latter will be applied to collisional 
decoherence in Section 5.2. This will allow us to evaluate the statistical weights of the 
pointer states in Section 5.3. 



5.1. Quantum trajectories and the orthogonal unraveling 

5.1.1. Quantum trajectories In the quantum trajectory approach, the wave function 
corresponding to a pure initial state \ip (0)) is propagated stochastically to generate 
pure state trajectories {\ipi (t))}, whose ensemble average recovers the solution of the 
master equation (fl"2"|) . i.e. 

exp(£t)\i>(0))(rJj(0)\=M(\rJ H (t))(i k (t)\) . (56) 

Such a stochastic process is called an unraveling of the master equation. A common 
unraveling is provided by the quantum Monte Carlo method [28, 29] which is based 
on a piecewise deterministic process. Here, one realization of a trajectory consists of 
smooth deterministic pieces generated by an effective (non-Hcrmitian) Hamiltonian, 
in our case 

ih f°° 

Hcff = H- y y dgLjL,, (57) 

which are interrupted by random jumps. The jumps occur with the rate 

r q = (L+L,) , (58) 

an expectation value with respect to (t)), and they are effected by the operators 
J 9 = L,/V?V (59) 

The quantum Monte Carlo method is not the only stochastic process satisfying (f5Gl) . 
In fact, there exists an infinite set of these stochastic processes, because the convex 
decomposition of the density matrix on the left hand side of (|56p is not unique. For 
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instance, other unravelings can be obtained from the quantum Monte Carlo method, 
by noting that the generator C does not uniquely fix the Lindblad operators L q and 
the Hamiltonian H [21,29]. This is due to the fact that the master equation (jT2J) is 
invariant under certain transformations of the Lindblad operators, such as the addition 
of a complex multiple z q of the identity I, 

\- q ^L' q = L q + z q \. (60) 

In the latter case, also the Hamiltonian must be transformed as 

1 f°° 

H^H' = H + -/ dq{z* q l q -z q ll) , (61) 
in order to assure the invariance of the master equation. 

5.1.2. The orthogonal unraveling To obtain a different, but again piecewise 
deterministic unraveling, we now make the choice z q = — (L„). It then follows that 
the deterministic pieces of a sample path are given by the solution of the nonlinear 
equation (|16p (or equivalently the corresponding projector equation (jlip - ). The jumps 
occur with the rate 

r g — C-gLqr) — (L^)(L g ) , (62) 

and are caused by the nonlinear operators 

J g = (l q -(L q ))/^. (63) 

As a distinctive feature, the states \ip q (t)) — J q \ijj (t)) into which the system may 
jump are orthogonal to the original state \ip(t)}, thus justifying its naming. (The 
states \ip q (t)} are not necessarily mutually orthogonal, though.) To our knowledge, 
this unraveling was first noted by Rigo and Gisin [16], although, it has not been studied 
numerically so far. 

A related unraveling, which is also referred to as the 'orthogonal unraveling', was 
introduced by Diosi [15,30]. Here, the deterministic pieces of the evolution are as well 
generated by (|16[) . However, the states \tp q (t)} into which the system may jump are 
obtained differently, as the eigenvectors of the Hermitian operator 

(1 - P^ (t) ) C (P m ) (1 - P^ (t) ) , (64) 

where P^,(t) = \ip (t))(ip (t) \. As a consequence, these states are also mutually 
orthogonal (in finite dimensional systems). Since the orthogonal unraveling of [15,30] 
requires the diagonalization of the operator (|64p . it is much more involved than the 
one defined by (fT6|) . (f62|) and (f63|) . which is why we will use the latter in the following. 

5.1.3. Pointer states and the orthogonal unraveling As mentioned above, all 
unravelings are equivalent in the same sense as the different convex decompositions 
of the density matrix e ct po. Note, however, that a preferred set of pure states 
- the pointer basis - may be singled out through the environmental coupling. In 
that case, those unravelings are distinguished which generate for any initial state an 
ensemble of these state independent projectors. An unraveling will do this job if (a) 
its deterministic part exhibits stable fixed points or solitons P a — |7r a )(7r Q |, which (b) 
are characterized by a vanishing jump rate, r (P Q ) = . In that case, the sample paths 
of the process will end up in one of the states P Q , by all means. Hence, the ensemble 
mean is of the form (J3J), such that the fixed points or solitons P Q can be identified 
with the pointer states of the system. 
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For the case of collisional decoherence the orthogonal unraveling fulfills the 
aforesaid conditions, (a) Its deterministic part is given by (TIT)) which exhibits the 
solitonic solutions ir (x) shown in Fig. 2. We shall see explicitly in Section 5.2.1 that 
these states are attractive fixed points, (b) We will show in Section 5.2.2 that the jump 
rate (|62[) vanishes for the solitons it (x), which finally demonstrates that our 'candidate' 
pointer states ir (x) are genuine pointer states P Q in the sense of (1). Moreover, this 
shows that the orthogonal unraveling is a very efficient numerical scheme for the long 
time solution of the master equation (TT3|) , since the state is no longer affected by the 
stochastic part, once it has turned into the soliton, and the trajectory is therefore 
more easy to integrate. 

We note that the orthogonal unraveling is not the only stochastic process which 
generates the ensemble of pointer states. In fact, there is a diffusive unraveling [31], 
which also involves (|llj) as its deterministic part. It was applied in [6] to investigate 
pointer states in a linear model. 



5.2. Unraveling collisional decoherence 

We now apply the described orthogonal unraveling to collisional decoherence (fT3|) . 
first evaluating the deterministic part (fl6|) of the stochastic process in Sect. 5.2.1 and 
then the stochastic one (|62p and ([63]) in Sect. 5.2.2. 



5.2.1. Deterministic evolution Applying (|16p to the case of collisional decoherence 
yields the soliton equation (fTT)) discussed in Sect. 3.2. We will now further simplify 
this equation, by considering initial states 

N 

V (x) = X>( )&M) , (65) 

2=1 

which are superpositions of non-overlapping wave functions fa (a;, 0), 

fa(x,0)4,*^(x,0) =0. (66) 
The latter are assumed to be localized in the sense that 

< < s- - (67) 

where c 2 . and a\, denote the variances of the distributions \fa (x, 0)| 2 and G(q), 
respectively. Under this assumption, which will be justified at the end of this section, 
one can extract a system of evolution equations for the time evolution of the coefficients 
in (E3, 



! I N N 

i=i j,k=l 



dt 



Here, the matrix = F (xi — Xj) is obtained from the localization rate (|15p. where 
the Xi = (x)0 4 denote the mean positions of the constituent wave functions fa (x,t). 
The latter evolve according to 

d t fa(x,t) = - ^—.dlfa{x,t) + fa{x,t)k\\faf\ (x,t) 
2mi L J 

N 

+ fa(x,t) l c i(*)| 2 7y(M) , (69) 
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where A is defined in (fl~8"|) and 7^ is a rate of the order of 7, 

7y (x, t) ee |^| 2 * F (x, t) - |0/ * F (a:, t) + Fy . (70) 

Let us now verify that ^f t {%) = $Zi=i c * CO 0* 0^ 0> with c i (i) an d & i) solutions 
of (155)) and (fBT))) . evolves according to (fT?]) . First, we note that the assumption (fBT)) 
leads to the approximation 

d^l^^V^-e^^, (71) 

for all q contributing appreciably to integrals weighted with the momentum transfer 
distribution G(q). This, in turn, implies 

Fjk = F (xj — Xk) (72) 
&x\^(x)\ 2 (|^| 2 *F) (x) . 

Hence, we find that dx \ <j>i (x)\ 2 (\4>i\ 2 * F^j (x) ~ 0, and one thus obtains 

a[|0 2 | 2 ] (x,t) = - (|<^| 2 *F) (x,t) . (73) 
Now, consider the time derivative of {x), which gives 

h N ( N \ 

i=l \i=l / 

N N \ 

£| Ci | 2 (|^| 2 *F)- £ |c/|c fc | 2 F,- fc , (74) 

where we dropped the arguments for brevity. In (|74|) . we used (|68l) - ([70|) . ([73]) . the fact 
that Fa — 0, and the normalization condition X^<=i l c »| 2 = 1- Now, we replace Fjk in 



74|) by the right hand side of (|72j) . and use (|66|) . which yields 

h 



d t ^t{x) = - —d^ t (x)-* t {x) (\V t \ 2 *F(x) 

- ^ dx|* t | 2 (y) (m 2 * F) (y)) . (75) 

J — OQ * 

Finally, by using (|15p . we obtain 

fit** (x) = - ^-d 2 x ^! t + 7 tf t (z) (76) 

x ^ t | 2 *G(x)-| °dy|^| 2 (y)(|* t | 2 *G) (y)) , 

which confirms that 'J't (x) evolves according to (fTT)) . 

Let us now discuss the solution of Eq. for the time evolution of the coefficients. 
We first consider situations where the wave packets <pi (x) are sufficiently far apart such 
that the localization rate is saturated, i.e. F^ = 7 (1 — Sij) . Under this assumption, 
reduces to the equation 



d ci(t) = -7(£M*)I 4 -M f )l 2 h(*)> (77) 



df 
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Figure 6. Numerical solution of H68H for N = 3. The rr-axis gives pi = |ci| 2 , the 
y-axis p2 = | C2 1 2 , and | C3 1 2 is fixed by normalization. Left: trajectories indicating 
the flow into the stable fixed points |cj| = Si >n ,(n = 1,2.3). Right: regions of 
attraction of the stable fixed points; the area denoted with n is the region of 
attraction of the fixed point |c;| = 5j in . 

which was already studied in [32] in the context of a discrete model for quantum 
measurement. It is shown there that all stable fixed points of (fT7| have the form 
\<h\ = <5;,n , and that the particular fixed point \ci\ = 8i tTn , with 

m = argmaxj (jc; (0)| 2 ) , (78) 

is approached monotonically, i.e. the component with the largest initial weight wins. 
This behavior is visualized in Fig. 6 which was obtained by solving (|TT[) numerically 
for the case N = 3. Here, the x- and the y-axis indicate the weights pi — |ci| and 
Vi = \ c 2\ 1 respectively. The plot on the left hand side shows various trajectories, 
illustrating in particular the fixed points. The plot on the right displays the regions of 
attraction of the stable fixed points |cj| = Si. n , in agreement with the criterion {78]). 
For instance, area 1 highlights the region of attraction of the fixed point \ci\ — 5i t \. 

Figure 7, on the other hand, depicts a scenario where the wave packets <f>i (x) are 
close together such that the localization rate is unsaturated, i.e. i<y ^ 7(1 — %). 
Here, we choose N = 3 wave packets with non-equidistant position expectations, 
(x±, X2, X3) crc/h = (1-4, 1.3, 0.8). Similarly to the saturated case, we observe that all 
stable fixed points of ([68]) have the form \ci \ = 5i. n . However, the regions of attraction 
are deformed such that criterion l|78j) is no longer valid, and the fixed points are not 
necessarily approached monotonically. 

To see that |cj| = 8i, n are stable fixed points of assume that the coefficients 
are close by, i.e. \cn (t)\ = 1 — e (t) with e (0) -C 1. It follows from ([68]) that 

i (*) = - E \°i W I' + ( £2 ) < ' ( 79 ) 

3 

and hence, \a (t — > oo)| — 8i >n . 

The knowledge of the fixed points of the coefficients allows one to discuss the 
asymptotic evolution of the initial state shown in Eq. (|65j) . Since the coefficients Cj 
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with j to tend to zero asymptotically, it follows that 

(t -f oo)) = \<p m (t -> oo)) , (80) 

for a specific to (which is given by |78|) in the saturated case) . The asymptotic behavior 
of the wave packets \fa n ) can, in turn, be predicted from Eq. (|69|) . Since the Cj-t m 
vanish for large times, the coupling term given by the last summand in (|69p vanishes 
as well, implying that the time evolution f|69|) of \fan) is asymptotically equal to the 
soliton equation (fT?)) . Therefore, in the absence of stochastic jumps, the state |\&o) 
evolves into that solitonic solution 7r m (x) of (fT7]) which is associated to the initial 
wave packet \(j> m (0)). 

It should be mentioned that Eqs. (|S5|) and (|69[) for the coefficients q and the 
constituent wave packets fa are not completely decoupled, since depends on the 
matrix i^j = F (xi — Xj ) which contains the position expectations Xi of the wave 
packets fa. However, the position expectation follows the classical trajectory for 
sufficiently large k's, such that can be solved without knowing the solution of 
dMJ). 

Let us now discuss the validity of the assumption of small position variance 
(|67p . and the ensuing approximation (|7ip . It can be justified by our observation in 
Section 4.3, that the dimensionless pointer width o-^cjq /ft is a function of the parameter 
k = GqI mJvy only, 

Ctt^ = -7— haioc, with a loc = 0.4 . (81) 

fi 4ai oc 

Thus, for all k <^ 4aj 2 oc ~ 1 we find that the position variance <j\ of a pointer state is 
one order of magnitude smaller than the reciprocal width of the momentum transfer 
distribution G (q), 

h 2 h 2 h 2 
°l * aL— - 0-2— < 2tt— . (82) 
a G a G a G 
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The above relation for the width of the pointer state is sufficient to justify the 
approximation J dxe lqx / h \ir (x)\ 2 ~ g^M^/^ ag we checked numerically, by using 
the solitonic solution of (6). The relative error is less than 2% for q 6 [ — 2(7(3, 2o~q] 
and k < 10~ 3 . 



5.2.2. Stochastic part Upon inserting the Lindblad operator L q = y/jG (q)e' lqx into 
(|63 jl we find that the jump operator takes the form 

J„ = Afq (e i9x/ri - {e iqx/h )^j , (83) 

with normalization Af q = (l - \{e iq */ A )\ 2 ) ' . We again consider states of the form 
(|65[) which are superpositions of non-overlapping (|66[) and localized (|67p wave packets 
(a;). Under this assumption, one can evaluate the expectation value in (|83|) 

N 

(e < «"/ A )*= 5^|c,-| a e <,x * /n , (84) 
»=i 

such that the state 'J'q (x) = J 9 'I' (x) into which the system may jump takes the form 

(N \ N 

eiqx /H _ J- | Cj | 2 e W> £ c, 0l (a; ) . (85) 
t=l / i=l 

Later we will choose the initial wave packets 4>i (x) to be solitons TTi (x). Let us 
therefore assume that the fa's form a basis, such that ^ q (x) can be represented as 
\& q (x) = Yli c i (?) i'i ( x ) ■ Then the transformed coefficients Ck (q) can be evaluated 
by the overlap Ck (q) = (<pk\-iq\^)- Using (|66[) and (|71[) this leads to the following 
expression for the redistribution of the coefficients due to an orthogonal jump 



c k (q) = M q [e iqx ^ h - ei9Xi/ ^ c * . 



(86) 

Similarly, one can evaluate the rate (|62|) associated to the jump operator of collisional 
decoherence, 

r q = 1 G{q)(l-\{e i ^ n )f) . (87) 
The above approximation (|84[) further simplifies this expression, 

N 



r q = lG{q) 1- \ c oV Vi.V<- 



j,k=l 

One observes that this rate vanishes for the stable fixed points |c,-| = Si >n , indicating 
that the quantum trajectories of the orthogonal unraveling evolve into the pointer 
states, the solitonic solutions of (fT"T)) . Moreover, we note that the original stochastic 
process (ITB1) . ([6"2j) . (f()3"j) - which is defined in the infinite dimensional Hilbert space 
of the system - has been reduced to a stochastic process in (D , demonstrating the 
efficiency of the orthogonal unraveling. However, due to the finite pointer width (|37[) 
the exact expression for the jump rate (|87[) does not vanish identically, although it is 
very small compared to 7. For instance, the numerically obtained soliton displays a 
strongly suppressed total jump rate r to t = / dqr q of r to t/l = 7 x 10~ 3 for n = 10~ 3 , 
while the superposition state decays with the rate r to t — 7- This implies that the 
solitons are not perfect pure state solutions of the master equation (|13p , though the 
loss of purity is small. 



Pointer basis induced by collisional decoherence 
1 1 <- 



22 



0.8 




0.4- 



0.2 



0.5 



1 t 1.5 



Figure 8. Quantum trajectory generated by 168t . J93D and (1941 with N = 2. 
The solid line depicts the evolution of |ci(t)| 2 , while the dashed line shows 
|c2 (t)| 2 - Since there is an odd number of jumps (three jumps in this example), 
the trajectory evolves into the fixed point |cj| = Si x. 



5. 3. The statistical weights of the pointer states 

The previous section showed that the orthogonal unraveling of an initial superposition 
state subject to collisional decoherence can be reduced to a stochastic process with 
respect to the corresponding coefficients. In particular, this applies to the case where 
the initial state is a superposition of pointer states, 

N 

|* >=5>M°)>- (89) 

i=l 

Thus we can now verify, by using the discrete process defined by Eqs. and 
([55]). that after decoherence the statistical weights of the pointer states are given by 
the overlap of the initial state with the initial pointer states. More specifically, this 
demonstrates that the initial state (x) evolves into the mixture 

N 

p(x,x') = (x)tt* (x') , (90) 

i=l 

where the statistical weights are given by the overlap 

Pl = |(*oMx,0))| 2 . (91) 

We first present an analytic proof of the above for N = 2 . The general case, N > 2 , 
is then treated numerically in the following section. 

5.3.1. Superposition of two localized states We consider the expectation value for the 
coefficients after a jump, that is 

/•OO 

(c k (q)) G := dqG(q)c k (q), k = 1,2. (92) 
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Upon inserting (|86|) one obtains (c\ {q))c — AP |c2 1 2 Ci , with the normalization con- 
stant AT' = (N q (e^/ h - e l ^' h )) G , see ®. Using |(ci (q)) G \ 2 + |(c 2 (<?)) G | 2 = 1, 
we find | J\P | = 1/ (|c 1 C2 1) which implies 

\(ci(q)) G \ = |ca| • (93) 

This shows that after an average jump the moduli of the coefficients are simply 
interchanged. This property (which does not hold for N > 2) makes the stochastic 
process analytically tractable, not least because the dynamics is independent of the 
phases of the coefficients. Since the deterministic part of the evolution is 

monotonic for N = 2, a trajectory starting from \c\ (0)| < 1/2 will end up in the 
state \ci (oo)| = (5,.i if and only if an odd number of jumps occurs in the process. This 
is demonstrated in Fig. 8. Crucially, the jump rate r to t (t), 

I* oo 

not (t) := / dqr q (t) 



= 2F(x 1 -x 2 )\c 1 (t)\ 2 \c 2 (t)\ 2 , (94) 

is unaffected by the jump (|93p at all times, since it is invariant under interchanging 
the coefficients. Hence, the time dependence of the jump rate is identical for all 
trajectories, which, in turn, means that the number of jumps follows an inhomogeneous 
Poisson process. Therefore, the probability for an odd number of jumps, which is equal 
to the statistical weight p\ of the pointer state m (x) , is given by 

Prob(odd) = (l-e- 2 ^°°V 2 ) , (95) 



with fi (t) — J Q * dr r to t (t) the integrated jump rate. The latter can easily be evaluated 
by noting that (|68|l can be written for N = 2 as 

2F(x 1 ~-x 2 )\c 1 (r)\ 2 \c 2 (r)\ 2 = \ A In (l - 2 \c x (r)| 2 ) . (96) 
Upon inserting this result into (|94p we obtain the integrated jump rate: 

^°°)=r dT ^ in ( i " 2ici(T)|2 ) 

= -ln(l-2| Cl (0)| 2 )/2. (97) 

Noting (f9"S")) we thus find the probability for an odd number of jumps 

Prob(odd) = |d (0)| 2 . (98) 

This finally confirms that the statistical weights of the pointer states are indeed given 
by the expected overlap (f9"Tj) . 



5.3.2. Superposition of N > 2 localized states The stochastic process is much 
more involved if the initial superposition consists of more than two pointer states. 
Our numerical implementation of the stochastic process defined by (|68p. (f5o| and 
(|88p is based on a Metropolis-Hastings algorithm to draw the momentum transfer 
q in accordance with the rate , with G (q) a Gaussian. Each of the generated 
trajectories ends asymptotically in one of the fixed points corresponding to a pointer 
state, and we thus obtain a numerical estimate of the statistical weights by means 
of the relative frequencies fk,l^k^N,of the asymptotic states. To confirm 
the expected probability distribution pk = \ck (0)| 2 , we evaluate the relative entropy 
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Figure 9. Relative entropy H {fk\Pk) of t ne numerically obtained distribution 
of pointer states f). with respect to the expected distribution pk = |cfc| 2 as a 
function of the number of trajectories n generated in the simulation. The plot 
indicates that the pointer states are distributed according to the initial overlap 

|c fc | 2 . 



H (fk\Pk) between these two distributions. Figure 9 shows the result for a random 
initial state with TV = 5 as a function of the number of trajectories n, indicating 
convergence to zero. In addition, we found for 100 random initial states, with random 
2 < N < 11, based on 10 4 trajectories that the relative entropy was always less than 
4 x 10~ 3 . This holds both for cases where the initial wave packets 7T; (x) are far 
apart such that the localization rate is saturated, F [xi — Xj) ~ 7, and for situations 
where the wave packets are located close together such that F (xi — Xj) < 7. We 
conclude that the asymptotic trajectories are indeed distributed according to the 
expected overlap (|9"Tj) . 

As an alternative confirmation of the statistical weights, we performed a % 2 -test. 
Similar to the treatment above, 100 random initial states {^l 1 ^ i ^ 100}, with 
random 2 < N < 11, were drawn by the simplex picking method. For each random 
state, n = 100 trajectories were generated, each of which ends asymptotically in one 
of the pointer states. Using the observed relative frequencies , 1 ^ k ^ N, of the 
pointer states, we evaluate 

2 W/*-M0)| 2 ) 2 

X 2 = nY± 5-^-, (99) 

ti M0)l 2 

for each random state. In order to verify that the pointer states are distributed 
according to | Ck (0) | 2 , the set { xf } must be shown to be sampled from a ^-distribution 
with N — 1 degrees of freedom. Comparing the set {xf} with the a-quantiles Q 
(denoted by Q a ) of the corresponding ^-distribution, a typical run shows ten cases 
where x\ > Q0.9 , one case where x\ > Q0.99, but not a single case where xf > Q0.999, 



X For instance, the 0.9 quantile is the value such that 90% of the samples lie below Qo.9- 
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as one expects if the {xf} are X 2 -distributed. Like above, this confirms statistically 
that the asymptotic trajectories are distributed according to the expected overlap 
(EH). 

6. Conclusion 

In this article, we related the nonlinear pure state equation discussed in [6-8] to a 
specific orthogonal unraveling of the collisional decoherence master equation. This 
puts into evidence that the dynamics of a particle in an ideal gas environment can 
be represented by an ensemble of pure state trajectories which evolve into spatially 
localized pointer states. For sufficiently strong collisions with the background gas, 
these solitonic wave packets move according to the classical equations of motion, thus 
explaining the emergence of classical dynamics within the quantum framework. 

Once the pointer state is reached by an individual quantum trajectory, the latter 
is no longer affected by the stochastic part of the unraveling, such that the integration 
of the trajectory is reduced to the solution of the classical equations of motion. This 
suggests that the orthogonal unraveling is an efficient algorithm for the long time 
solution of master equations which exhibit a pointer basis. On the other hand, also 
the short time solution turns out to be efficient, since the orthogonal unraveling can 
be reduced (under appropriate assumptions) from an infinite dimensional unraveling 
to a stochastic process in (D . 

Future studies might consider the emergence and dynamics of pointer states in 
dissipative quantum systems. We note that the present work relies on the model 
of pure collisional decoherence which does not describe long time effects such as 
dissipation or thermalization. It would certainly be worth to determine the pointer 
states of a more involved model such as the quantum linear Boltzmann equation 
presented in [33-35]. For large mass ratios between the test particle and the gas, one 
expects that the pointer states then evolve according to a Langevin equation, thus 
explaining the emergence of classical Brownian motion within the quantum framework. 

We thank B. Vacchini for helpful discussions. The work was supported by the 
DFG Emmy Noether program. 
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